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geologic unit [3], A uniform, globally higher heat flow should not 
produce such localized sapping valleys. It is for these reasons that 
we invoke localized hydrothermal systems for martian valley gen- 
esis. Such hydrothermal systems would be localized in surface 
extent, yet (as the simulations suggest) draw groundwater from 
great distances while focusing outflow into relatively small regions. 
Th e hydrolherm a 1 d i sch a rge wou Id a 1 so pre fe rent i a 1 1 y produ ce 1 a nd - 
forms associated with groundwater outflow — the sapping valleys. 
However, snowfall melting in hydrothermal areas and sublimating 
elsewhere might also play a role in producing such a distribution [3] . 

The principle Finding of our numerical model is that magmatic 
intrusions of several 10 2 km 3 provide sufficient volumes of ground- 
water outflow over the timescales (several 10 5 yr or more) needed 
to form fluvial valleys [8]. Calculated discharges are robust. Subsur- 
face inhomogeneilies, local impermeable caps, and uncertainties in 
porosity all affect the discharges at the 20% level or less. The param- 
eter with the single greatest effect on the calculated discharges is the 
subsurface permeability. Permeabilities between lOand lOOOdarcy 
provide sufficient quantities of groundwater out flow to form fluvial 
valleys. Lower permeabilities require larger intrusion volumes to 
produce the same discharge. However, in the range of permeabi lilies 
expected for basaltic rock, there is no difficulty in producing signifi- 
cant groundwater outflow. 

We have estimated the eroded ^volumes of two of the best- 
developed valley networks (Parana and Warrego Valles) in the 
heavily cratered terrains. These values (H.C.T) are compared with 
our estimated valley volumes on martian volcanos in Fig. 1. Esti- 
mates of martian valley erosion can be combined with terrestrial 
fluvial erosion rates to obtain estimates of the total volume of water 
required to form each set of martian valleys. Some ratios of water 
volume to eroded volume for Mars are as lpw as 2 or 3 to 1 [9]. 
However, based upon our own study of fluvial erosion on volcanic 
landscapes, we find ratios as large as 1000: 1. The total water volume 
using each ratio is shown for each valley group in Fig. 2. For each 
locality the lower bar represents the uncertainty due to valley side- 
wall slopes while using a water-to-eroded-volume ratio of 3: 1, the 
upper barusing a ratio of 1000: 1 .Horizontal lines in Fig. 2 illustrate 
the cumulative discharge of hydrothermal systems associated with 
50-,500-,and 5000-km 3 igneous intrusions. Therefore we conclude 
that hydrothermal systems can provide the volumes of groundwater 
outflow needed to form martian valley networks and can provide an 
alternative to rainfall from a warm, wet early Mars. 


There are several lines of evidence that suggest early Mars was 
warmer and welter than it is at present [1]. Perhaps the most 
convincing of these are the valley networks and degraded craters 
that characterize much of the ancient terrains. In both cases, fluvial 
activity associated with liquid water is believed to be involved. 
Thus, Mars appears tohavehad a warmer climate early in its history 
than it does today. How much wanner is not clear, but a common 
perception has been that global mean surface temperatures must 
have been near freezing — almost 55 K warmer than at present. 

The most plausible way to increase surface temperatures is 
through the greenhouse effect, and the most plausible greenhouse 
gas is C0 2 . Pollack et al. [2] estimate that in the presence of the faint 
young Sun, the early martian atmosphere would have to contain 
almost 5 bar of C0 2 to raise the mean surface temperature up to the 
freezing level; only 1 bar would be required if the fluvial features 
were formed near the equator at perihelion al maximum eccentric- 
ity. However, these calculations now appear to be wrong since 
Kasting [3] has shown that C0 2 will condense in the atmosphere at 
these pressures and that this greatly reduces the greenhouse effect 
of a pure C0 2 atmosphere. He suggested that alternative greenhouse 
gases, such as CH 4 or NH 3 , are required. 

In lit is paper, we approach the early Mars dilemma from a 
slightly different point of view. In particular, we have constructed 
a model for the evolution of C0 2 on Mars that draws upon published 
processes that affect such evolution. Thus, die model accounts for 
the variation of solar luminosity with time, the greenhouse effect, 
regolith uptake, polar cap formation, escape, and weatheringlWe 
initialize the model 3.8 G.y. ago with a specified C0 2 inventory and 
iheiLmarch it forward in time to the present epoch. The model 
partitions C0 2 between its various reservoirs (atmosphere, caps, 
Tegolith, carbonates, and space) according to the thermal environ- 
ment predicted by a modified version of die GTerasch and Toon [4] 
energy balance climate model. The goal is to determine if it is 
possible to find an evolutionary scenario dial is consistent with early 
fluvial activity, and that arrives al the present epoch with the initial 
C0 2 partitioned intoits various reservoirs in plausible amounts. Our 
early fluvial activity criterion is that global mean temperatures must 
be at least 240 K at flie beginning of the simulation^ our current 
reservoir criteria is that the atmosphere must hold about 7.mbar of 
C0 2 , the caps several millibars, and tire regolidi 300 mbar. Wg do 
not constrain the Final size of the rock reservoir. ~ _ 

WeJTInd no evolutionary scenario that salaries these criterion 
_ whenusing published estimates of the processes involved. The main 
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Fig. 1. Temperature (top) and pressure evolution (bottom) for total initial 
abundance of 1.5 bar. The temperatures shown are for the surface at equator 
(solid line with circles), pole (solid line with squares), global mean (dashed line 
with triangles), and frost point (dotted). The C0 : reservoir pressures are 
atmosphere (solid line with squares), caps (dashed line with circles), regolith 
(dotted line with triangles), and carbonate (chain dotted line with plus sign). 


difficulty is making Mars warm early on. As Kasling has pointed 
out, only a stronger greenhouse or a brighter early Sun can help in 
this regard. However, we have found that if the greenhouse were 
stronger or if the Sun were brighter, then massive permanent caps 
would form as a result of a collapse in the climate system sometime 
between 1 .5 and 3,5 b.y. ago. An example of this collapse and the 
thermal history associated with it is shown in Fig. 1 . 

The collapse is a result of an unstable feedback between the 
poleward transport of heat by die atmosphere, the greenhouseeffect, 
and surface pressure. As surface pressure falls, heat transport and 
the greenhouse effect are reduced, die polar caps cool, surface 
pressure falls further, and so on. Gierascb and Toon [4] discuss this 
instabilhy in detail. In our model, die instability is set off by 
weadiering that removes C0 3 from the atmosphere at a rate that is 
exponentially proportional to temperature. Thus, die warmer early 
Mars is. die more likely collapse will occur. As much as 600 nibar 
goes into the caps when collapse occurs with the C0 2 coming from 
the atmosphere and die regolith. More importantly, at least 300 nibar 
survives to the present epoch — much more Uian appears to reside in 
the south residual cap. 


Collapse can be avoided if die polar albedo is significantly lower 
than the value we have assumed (0.75), or if the actual poleward 
heat flux is greater than that given by our simple parameterization. 
However, ip either case, the implication is that if global mean 
surface temperatures were at or above 240 K on early Mars, then a 
minimum total inventory of 2 bar of C0 2 is required, and at least 
70% of it has been sequestered as carbonate in near- surface mate- 
rials. On the other hand, if the fluvial features in the ancient terrains 
do not require global mean temperatures near 240 K and can be 
explained by phenomena that are not climate related, such as an 
elevated geothermal heal flux, then our model favors an initial C0 2 
inventory near 600 mbar. Of this initial C0 2 , most has gone into the 
regolith (300 mbar), modest amounts., into carbonates (130 nibar), 
even smaller amounts into the atmosphere (7 mbar) and caps 
(3 mbar), with the remainder having escaped, into space ( 1 60 mbar). 
Thus, it is crucial that we obtain belter constraints on the thermal 
regime required to form the fluvial features on early Mars, 
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We examine the effects of loss of Mars atmospheric constituents 
by solar-wind-induced sputtering and by photochemical escape 
during die last 3.8 b.y. Sputtering is capable of efficiently removing 
all species from the upper atmosphere, including the light noble 
gases; N also is removed by photochemical processes. Due to the 
diffusive separation by mass above the homopause, removal from 
the top of die atmosphere will fractionate the isotopes of each 
species, with the Iightermass being preferentially lost. For C and O, 
this allows us to determine the size of nonatmospheric reservoirs 
that mix with the atmosphere; these reservoirs can be accounted for 
by exchange with C0 2 adsorbed in the regolith and with H 2 0 in the 
polar ice deposits. We have constructed both simple analytical 
models and time-dependent models of the loss of volatiles from and 
su ppl y to the mart ian atmosphere . B o tli Ar and Ne requ ire cont i nu ed 
replenishment from outgassing over geologic time. 

For Ar, sputtering loss then explains the fractionation of 36 Ar/ 
38 Ar without requiring a distinct epoch of hydrodynamic escape 
(although fractionation of Xe isotopes still requires a very early 
hydrodynamic escape). For Ne, the currejil ratio of :2 Ne/ 20 Ne rep- 
resents a balance between loss to space and continued resupply from 
the interior; the similarity of the ratio to the terrestrial value is 
coincidental. For N, the loss by belli sputtering and photochemical 



